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Abstract — Q In this paper, collaborative use of relays to form 
a beamforming system with the aid of perfect channel state 
information (CSI) and to provide communication in physical- 
layer security between a transmitter and two receivers is in- 
vestigated. In particular, we describe decode-and-forward based 
null space beamforming schemes and optimize the relay weights 
jointly to obtain the largest secrecy rate region. Furthermore, the 
optimality of the proposed schemes is investigated by comparing 
them with the outer bound secrecy rate region. 

I. INTRODUCTION 

The open nature of wireless communications allows for the 
signals to be received by all users within the communication 
range. Thus, secure transmission of confidential messages is 
a critical issue in wireless communications. This problem 
was first studied in JT| where Wyner identified the rate- 
equivocation region and established the secrecy capacity of the 
discrete memoryless wiretap channel in which eavesdropper's 
channel is a degraded version of the main channel. Later, 
Wyner's result was extended to the Gaussian channel in |[3) 
and recently to fading channels in |4|. In addition to the single 
antenna case, secure transmission in multi-antenna models is 
addressed in J5) - (6). For multi-user channels, Liu et al. [7| 
presented inner and outer bounds on secrecy capacity regions 
■ for broadcast and interference channels. The secrecy capacity 
of multi-antenna broadcasting channel is obtained in lf8~l . 
Moreover, it's well known that that users can cooperate to form 
a distributed multi-antenna system by relaying. Cooperative 
relaying with secrecy constraints was recently discussed in 

GoJ-IEQ . 

In this paper, we study the relay-aided secure broadcasting 
scenario. We assume that the source has two independent 
messages, each of which is intended for one of the receivers 
but needs to be kept asymptotically perfectly secret from the 
other. This is achieved via relay node cooperation in decode 
and forward fashion to produce virtual beam points to two 
receivers. The problem is formulated as a problem of designing 
the relay node weights in order to maximize the secrecy rate 
for both receivers for a fixed total relay power. We assume 
that the global channel state information (CSI) is available for 
weight design. Due to the difficulty of the general optimization 
problem, we propose null space beamforming transmission 

' This work was supported by the National Science Foundation under Grant 
CCF - 0546384 (CAREER). 
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schemes and compare their performance with the outer bound 
secrecy rate region. 

II. Channel 

We consider a communication channel with a source S, 
two destination nodes D and E, and Al relays {i? m }^ =1 as 
depicted in Figure Q] We assume that there is no direct link 
between S and D, and S and E. We also assume that relays 
work synchronously and multiply the signals to be transmitted 
by complex weights to produce virtual beam points to D and 
E. We denote the channel fading coefficient between S and 
R m as g m € C , the channel fading coefficient between R m 
and D as h m € C, and the channel coefficient between R m 
and E as z m € C. In this model, the source S tries to transmit 
confidential messages to D and E with the help of the relays 
. It is obvious that our channel is a two hop relay network. 
In the first hop, the source S transmits x s which contains the 
confidential messages intended for both D and E to the relays 
with power i?[|x s | 2 ] = P s . The received signal at relay R m 
is given by 

Vv.ra — 9m%s ~t~ (1) 

where rj m is the background noise that has a Gaussian distri- 
bution with zero mean and a variance of N m , 

In the first hop, the secrecy rates for destination D and E 



lie in the following triangle region. 
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where Rd and i? e denote the secrecy rates for destination D 
and E, respectively. 

III. Relay Beamforming 
We consider the scenario in which relays are much more 
closer to the source than the destinations, and hence, the 
first-hop rate does not become a bottleneck of the whole 
system. Due to this assumption, we in the following focus on 
characterizing the secrecy rate region of the second-hop. We 
consider the decode-and-forward relaying protocol in which 
each relay R m first decodes the message x s ,and subsequently 
scales the decoded messages to obtain x r = w m Xd + u m x e , 
where w m and u m are the weight values. Xd and x e are 
independent, zero-mean, unit-variance Gaussian signals which 
include the confidential messages to D and E, respectively. 
Under these assumptions, the output power of relay R m is 

E[\x r \ 2 ] = E[\w m x d + u m x e \ 2 ] = \w rn \ 2 + \u m \ 2 (4) 

The received signals at the destination nodes D and E are 
the superpositions of the signals transmitted from the relays. 
These signals can be expressed, respectively, as 
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where n$ and n\ are the Gaussian background noise com- 
ponents at D and E, respectively, with zero mean and 
variance No. Additionally, we have above defined h = 
[h$,....h* M ] T ,z = [z{,....z* M ] T , w = K,...ui M ] T , and 
u = [ui, ...um} T ■ In these notations, while superscript * 
denotes the conjugate operation, ( ) T and (-)T denote the 
transpose and conjugate transpose , respectively, of a matrix or 
vector. From the transmitting and receiving relationship in (O 
and ©, we can see that the channel we consider can be treated 
as an interference channel with secrecy constraints studied in 
0. The achievable secrecy rate region is shown to be 
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In this paper, we address the joint optimization {w m } and 
{u m } with the aid of perfect CSI, and hence identify the 
optimal collaborative relay beamforming (CRB) direction that 
maximizes the secrecy rate region given by (0 and (0. Since 
the optimization problem above is in general intractable, we 
investigate suboptimal schemes. 



A. Single Null Space Beamforming 

In this scheme, we choose one user's (e.g., E) beamforming 
vector (e.g., u) to lie in the null space of the other user's 
channel. With this assumption, we eliminate the user _E's 
interference on D and hence D's capability of eavesdropping 
on E. Mathematically, this is equivalent to | Em=i h m u m \ 2 = 
W\i = 0, which means u is in the null space of h^. 

We further assume a fraction of total relay transmitting 
power P r is used for sending confidential message to D. Under 
these assumptions, we can solve the optimization problem in 
(0. The maximum Rd can be computed as 



i?d, m (h, z,P r ,a) 
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Here, we use the fact that (fL2l is the Rayleigh quotient 
problem, and its maximum value is as given in ( fT3b where 
A m aa;(A, B) is the largest generalized eigenvalue of the matrix 
pair (A, B). Note that we will also use X max (-) to denote 
largest eigenvalue of the matrix in later discussion. The 
optimum beamforming weights w is 
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where ip w is the eigenvector that corresponds to X max (NoI + 
aP r hhT,iVoI + aP r zz*) and <; is chosen to ensure 
™l P t™opt = aP r . 

Now we turn our attention to the maximization of R„ when 



w = w op4 . Note that N + \ J2™=1 
denoted by A^, Due to the null space constraint, we can write 
u = H^v, where denotes the projection matrix onto 
the null space of h^. Specifically, the columns of are 
orthonormal vectors which form the basis of the null space of 
ht. In our case, is an M x (M — 1) matrix. The power 
constraint u^u = vtH^H^v = < (1 — a)P r . 
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The maximum R e under this condition can be computed as 
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The optimum beamforming vector u is 
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where qi is a constant introduced to satisfy the power con- 
straint. Hence, secrecy rate region R s .b achieved with this 
strategy is 



< R d < R d , m (h,z,P r ,a) 
< Re < R e , m (h,z,P r ,a) 
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Note that we can switch the role of D and E, and choose 
w to be in the null space of zK In general, the union of region 
described in (|2T1 > and its switched counterpart is the secrecy 
rate region of single null space beamforming strategy. 

B. Double Null Space Beamforming 

In this scheme, we simultaneously choose the beamforming 
vectors for D and E to he in the null space of each other's 
channel vector. That is | Y2m=i h m Um\ 2 — h^u = 0, and 



Sm=i z rn.Wm\ 2 = z^w = 0. In this case, the channel reduces 
to two parallel channels. Since interference is completely 
eliminated, the secrecy constraint is automatically satisfied. 
Coding for secrecy is not needed at the relays. The channel 
input-output relations are 
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Now, we only need to solve the following problems: 
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Similarly as in Section IIII-AI we can easily find the secrecy 
rate region M.a,b for double null space beamforming as 
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where denote the projection matrix onto the null space of 
and is defined similarly as H^. 

C. TDMA 

For comparison, we consider in the second-hop that the 
relay only transmits secret information to one user at a time 
and treat the other user as the eavesdropper. We assume that 
relay uses a fraction of time to transmit Xd where (1 — a) 
fraction of the time is used to transmit x e . The channel 
now is the standard gaussian wiretap channel instead of an 
interference channel. It can be easily shown that the rate region 
Kjdmo is 

<R d < a log X max (N I + P r hh\ N I + P r zz^) (28) 

<R e < (1 - a) log A ma ,(^ I + ^zzt, N I + P r hh^) 

(29) 

IV. Optimality 

In this section, we investigate the optimality of our proposed 
null space beamforming techniques. Although the optimal 
values of w and u that maximize the rate region (0 and <[8J is 
unknown, we can easily see that the following rate region is 
an outer bound region of our original achievable secrecy rate 
region. 
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Again, this rate region should be maximized with all possible 



w and u satisfying ||w| 



lull 2 < P r . From the above 



expressions, we can see that this outer bound can be inter- 
preted as two simultaneously transmitting wire-tap channels. 
Fortunately, the optimization problem in this case can be 
solved analytically. With the same assumptions as before that 
||w|| 2 = aP r , ||u|| 2 = (1 — a)P r , we can easily show that 
the outer bound secrecy rate region M, 
relay beamforming system is 
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The expression for Rd and R e here coincide with the secrecy 
capacity of Gaussian MISO wiretap channel with 
transmit power levels aP and (1 — a)P, 

A. Optimality in the High-SNR Regime 

In this section, we show that the outer bound region M. outer 
converges to the proposed null space beamforming regions at 
high SNR. For the single null space beamforming scheme, 
the maximum Rd in ( TT3l ) has the same express as in d34l . 
and thus it is automatically optimal. R e in single null space 
beamforming has basically the same expression as that of R e 
in double null space beamforming with N replaced by N t . 
This difference is negligible as P goes infinity. Hence, we 
focus on double null space beamforming and show that in 
the high-SNR regime, the M. outer coincide with the double 
null space region described by d26l i and d27l >. In the following 
analysis, for simplicity and without loss of generality, we 
assume No = 1. From the Corollary 4 in Chapter 4 of 151, we 
can see that 

lim — A ma x(I + -Prhh t ,I + P r zz t ) = max|h + ^| 2 (36) 

P r — foo P r S 

where ifi is a unit vector on the null space of zt. Similarly, 
we can define ipx as a unit vector on the null space of ht. 
Combining this result with ([34-b and d35l l. we can express the 
region M. ou t er at high SNRs as 

< R d < log(aP r ) + log(max |h f ^| 2 ) + o(l) (37) 

0<R e < log((l - a)P r ) +log(max|z t VJi| 2 ) + o(l) (38) 

where o(l) — > as P r —> oo. On the other hand, double null 
space beamforming region satisfies 

/ M \ 
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Above, d40b follows from the observation that 
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where ip is a unit vector and is in the null space of because 
w is in the null space of z^. (l42l follows similarly. Thus, the 
outer bound secrecy rate region converges to the double null 
space beamforming region in the high-SNR regime, showing 
that the null space beamforming strategies are optimal in this 
regime. 



B. Optimality of TDMA in the Low-SNR Regime 

In this section, we consider the limit P r — > 0. In the follow- 
ing steps, the order notation o(P r ) means that o(P r )/P r 
as P r -> 0. 

A mox (I + P r hht,I + P r zzt) (46) 
= A mQX ((I + P r zzt)- 1 (I + P r hht)) (47) 
= \ max ((I - Przzt + o(P r ))(I + P r hht)) (48) 

- Xmax ((I - PrZ t z)(I + P r hh r )) + o{P r ) (49) 

= X max (I + P r (hht - zzt )) +o(P r ) (50) 
= 1 + P r A mox (hh t - zz^ + ( Pr ) (51) 

Combining this low-SNR approximation with ( l34l l and ( l35l >. 
we can see that the M. outer at low SNRs is 

< R d < logA^^ + aP.hht.I + aP.zz 1 ") 

= aP r A maa: (hh t -zz t ) + o(P r ) (52) 

< Re < \og\ max (I + (1 - a)P r zz t , I + (1 - aPOW) 

= (1 -a)P r A mM (zz t -hh f ) + o(P r ) (53) 

Note that d52l l and (|53T > are also the low-SNR approximations 
for the TDMA approach. Thus, the TDMA scheme can achieve 
the optimal rate region in the low-SNR regime. For the 
completeness, we give the lower SNR approximations for 
single and double null space beamforming as well. For single 
null space beamforming scheme, the low-SNR approximation 
of (EB is 

< R d < aP r \ max (htf - zzt) + o(Pr ) (54) 
< R e < (1 - a)P r /N t z^H^z + o{P r ) (55) 

while for the double null space beamforming scheme, low- 
SNR approximations of (l26l l and (|27T i are 

< R d < aP-htH^H^h + o(P r ) (56) 
0< R e < (l-a)Pz t H^H^ t z + o(P r ) (57) 

C. Optimality when the Number of Relays is Large 
It is easy to show that 

A ma x(I + aP-hht,I + aP r zzt) < X max (I + aP r hht) 

= l + aP r hth (58) 

Now, consider the function 

1 + aP-htH^H^h (59) 

which is inside the log function in the double null space 
beamforming Rd boundary rate d26i l. In our numerical results, 
we observe that when M is large and h and z are Gaussian dis- 
tributed (Rayleigh fading environment), (l58l and ( |59l converge 
to the same value. Similar results are also noted when R e in 
(|27| | is considered. These numerical observations indicate the 
optimality of null space beamforming strategies in the regime 
in which the number of relays, M, is large. 



4 



2.25 
2.2 
2.15 



R (bits) 



Fig. 2. Second-hop secrecy rate region <r^ = 2, a z = 2, P r = 1, M 
Lower figure provides a zoomed version. 



V. SIMULATION RESULTS 

In our simulations, we assume N m — No = 1, and {g m }, 
{h m }, {zm} are complex, circularly symmetric Gaussian ran- 
dom variables with zero mean and variances a 2 g , a\, and <j\ 
respectively. 

In Figures [2] and [3] we plot the second-hop secrecy rate 
region of different schemes in which we see M. outer D R s j, D 
Kd,6 ^ Ktdma- We notice that our proposed suboptimal 
beamforming region is very close to outer bound secrecy 
region M. outer - Furthermore, the larger the M, the smaller 
the rate gap between M. ou t er and our proposed null space 
beamforming schemes. Also, we note that increasing the 
number of relays, M, enlargens the rate region. Moreover, 
we can see that M = 15 is sufficient for the null space 
beamforming schemes to coincide with the M. outer - 

Next, we examine the null space beamforming's optimality 
in the high-SNR regime in Fig. [4] In this simulation, we 
can see that when the relay power is large enough, M. oute r 
coincides with the regions of our proposed null space beam- 
forming schemes as expected even M is very small. Finally, 
in Fig. [5] where relay power small, we observe that R oute r 
coincides with the rate region of the TDMA transmission 
scheme. Also, we note that the double null space beamforming 
has better performance than single null space beamforming at 
some operation points. This is mainly because N t is no longer 
negligible at very low SNR values. 



Fig. 3. Second-hop secrecy rate region cr^ = 2, cr z = 2, P r = 1, M = 15 




Fig. 4. Second hop secrecy rate region ct/j = 2, ct z = 2, P r = 100, M - 



Fig. 5. 
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Second hop secrecy rate region = 2, <r z = 2, P r = 0.001, M 



5 



VI. Conclusion 

In this paper, we have considered a DF-based collaborative 
relay beamforming protocol to achieve secure broadcasting 
to two users. As the general optimization of relay weights 
is a difficult task, we have proposed single and double null 
space beamforming schemes. We have compared the rate 
regions of these two schemes and the TDMA scheme with 
the outer bound secrecy rate region of the original the relay 
beamforming system. We have analytically shown that null 
space beamforming schemes are optimal in the high-SNR 
regime, and TDMA scheme is optimal in the low-SNR regime. 
In our numerical results, we have seen that our proposed 
null space beamforming schemes perform in general very 
close to outer bound secrecy rate region. We have numerically 
shown that when the number of relays is large, the null space 
beamforming schemes are optimal. 
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